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We study static, spherically symmetric black hole solutions of the Einstein equa- 
tions with a positive cosmological constant and a conformally coupled self interacting 
\ scalar field. Exact solutions for this model found by Martinez, Troncoso, and Zanelli, 

CN | (MTZ), were subsequently shown to be unstable under linear gravitational pertur- 

bations, with modes that diverge arbitrarily fast. We find that the moduli space of 
static, spherically symmetric solutions that have a regular horizon -and satisfy the 
weak and dominant energy conditions outside the horizon- is a singular subset of 
a two dimensional space parameterized by the horizon radius and the value of the 
scalar field at the horizon. The singularity of this space of solutions provides an 
explanation for the instability of the MTZ spacetimes, and leads to the conclusion 
J> ■ that, if we include stability as a criterion, there are no physically acceptable black 

hole solutions for this system that contain a cosmological horizon in the exterior of 
^ — , \ its event horizon. 
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I. INTRODUCTION 



When one considers possible fields interacting with a black hole, the simplest source of 
matter that one could naively take into account corresponds to a single real scalar field. 
However, when this field is minimally coupled and the spacetime is asymptotically flat, the 
so called no-hair conjecture [l|, H, [1] indicates that this class of black hole does not exist. 
Much effort have been focused on this problem and recent works dealing with this issue can 
be found in pj. Nonetheless, this conjecture can be circumvented in different ways as we 
show below. 

A black hole solution, where the scalar field is conformally coupled, i. e. when the 
corresponding stress-energy tensor is traceless, was found in [5|. In this three-dimensional 
black hole, the scalar field is regular everywhere and the spacetime is asymptotically anti-de 
Sitter because a negative cosmological constant is included. This black hole solution can be 



*Elcctronic address: gdotti@famaf.unc.edu. ar 
^Electronic address: gleiser@fis.uncor.edu 
^Electronic address: martinez@cecs.cl 



2 



extended by considering a conformal self-interacting potential. This was done in [6(, where 
exact black hole solutions are found for a minimally coupled scalar field and a one parameter 
family of potentials. A previous four-dimensional and asymptotically flat black hole [7J was 
reported back in the 70 's, but the scalar field diverges at the horizon. The presence of a 
cosmological constant allows to find exact four-dimensional black hole solutions, where the 
scalar field is regular on and outside the event horizon (si, 0, [Io|, 11]. Numerical black hole 
solutions can also be found in four (H, El Q E, and five dimensions [17f. Further 



exact solutions in the context of low energy string theory were found in [18] . 

Some interesting aspects of these black hole solutions are studied in [l9[]. In particular, 
the analysis of stability against linear perturbations for the de Sitter conformally dressed 
black hole jl] done in 0], is relevant for the discussion presented here. 



In this work, we study the space of static, spherically symmetric solutions of the Einstein 
equations with a positive cosmological constant and a conformally coupled self-interacting 
scalar field (MTZ model). Conformally coupled scalar fields in General Relativity have 
been used to model quantum effects in semiclassical theories [2l[. This model has a 
well posed initial value formulation [22!], and was shown to reproduce better -than the 
minimally coupled scalar field- the local propagation properties of Klein Gordon fields 
on Minkowski spacetime [231 ] . Our interest, however, comes from the fact that this 
model allows non trivial static black holes solutions [8|. These solutions belong to a 
restricted class (eq. (|T2|) below) of spherically symmetric static spacetimes, and are given 
in equations (ITS]) (ffll) (solution MTZ1) and (Tf5]) (solution MTZ2). Note that a generic 
spherically symmetric static spacetime metric admits the local form (110]) . In this work we 
address the following question: Are there other static, spherically symmetric black hole 
solutions for the MTZ model, satisfying the dominant and strong energy condition be- 
tween the event and cosmological horizon, besides MTZ1 and MTZ2? Using a combination 
of analytical and numerical methods we conclude that the answer to this question is negative. 

The paper is organized as follows: in the next section we review the MTZ model and 
skecth the derivation of the MTZ1 and MTZ2 solutions. We also prove that MTZ1 is 
unstable under spherically symmetric gravitational perturbations. MTZ2 had already been 
found to be unstable under spherically symmetric gravitational perturbations, this being our 
original motivation to study the space of spherically symmetric static solutions of the MTZ 
system. This is done in Section lllfl where the full set of Einstein and scalar field equations 
is reduced to a second order ODE system. In Section [TV] we analyze the restrictions that 
the existence of a regular event horizon impose on the solutions, if we also require that the 
energy-momentum tensor satisfies appropriate energy conditions. Acceptable local solutions 
are found to be parameterizable with the horizon radius r and the value of the scalar field 
at the horizon, a := 0(r o ). The subset of allowed values is displayed in Fig. [2J To address 
the issue of the global behavior of these local solutions, the field equations were numerically 
integrated away from the horizon. Some illustrative examples are presented in Section |V] 
where the different behaviors as we move away from the event horizon are shown. We find 
that solutions that satisfy the energy conditions near the event horizon contain, in general, a 
coordinate singularity for some finite r outside the event horizon. We show that the isotropy 
spheres reach a maximum radius r at this point and contract as the proper distance from the 
horizon further increases. This explains why r is not an appropriate coordinate in this region. 
We provide an appropriate coordinate extension in Section IV Al A numerical integration 
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beyond the coordinate singularity, described in Section IV Bl suggests that, generically, the 
metrics contain a curvature singularity (the energy density diverges) at some finite proper 
distance from the extension point. 

Mostly for completeness we include in Section IV CI an analysis of the metrics that violate 
the energy conditions. A compilation of the main results, together with some final comments 
and our conclusions are given in Section [VI] 



II. THE MTZ MODEL 



In the MTZ model [8[ the action for gravity conformally coupled to a scalar field with 
a quartic self-interaction potential and an electromagnetic field F^ v is given by, 
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where a is a coupling constant. Variation of this action with respect to the metric, scalar 
field and Maxwell potential gives the following set of Euler Lagrange equations: 
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where the stress energy tensors are 



\g,u9 a(S d a( pd^ + ~ [g, u D - V M V, + G, u ] 2 - ag^ A (2d) 
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(2e) 



Under conformal transformations — > Q" 1 ^, F^ —>■ F^, g^ v — > fi 2 ^, equations pb|) - fl2c|) 
are invariant, and T^ u — > Q~ 2 T[^ l/ , T™ — > Q~ 2 T^ U M . This is the motivation behind the 
choice of the non minimal coupling and quartic self interaction of the scalar field. 
Note that the trace of T$ vanishes on shell: 
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(3) 



whereas T EM := Tzf 1 g^ u vanishes identically. Thus, taking the trace of Eq. (|2al) gives 



R = 4A 



(4) 



We should stress here that does not follow from l[2a]) alone, but from the system l[2a]) - 
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It is interesting to comment on those solutions of the field equations (|2J) for which = O , 
O 7^ a constant (we are not interested in the pure Einstein-Maxwell case = 0). In this 
case, the system (|2D reduces to: 

(l - ^pj G, v + (A + a^)g, u = T% M (5a) 

R + 2Aa<p Q 2 = (5b) 
Vi^ = 0. (5c) 



Taking the trace of (j5aj) and using fl5bl) gives 



(9a) 

(9b) 
0, (9c) 



Thus ( |5al) takes a simple form: 



(6) 



(l + 3^) [GW + A^l = ^ (7) 

and (|5bp gives again (T4]). Note that these are Einstein-Maxwell equations with an effective 
Newton's constant G e ff = (l + G [8], thus the case where (l + t^-) < (negative 

G e //) is somewhat pathological because it is equivalent to having repulsive gravitational 
forces j8|. 

The theory with a coupling constant a tuned with the cosmological constant as 

A 

is particularly interesting, since it seems to admit a wider set of solutions. We will call these 
theories special from now on, and call a ^ —A/36 theories generic. For special theories and 
constant scalar field configurations, O 2 = 6 and the field equations (jHJ) become 



= 



R-AA 

v ± av 



Note that the Euler-Lagrange equation for the metric, eq. (|9a|) . gives no information 
about the metric, but implies F^ u = 0. This does not mean that the gravitational field is 
unconstrained, as one might first be lead to think, since the Euler-Lagrange equation for 
the scalar field forces R = 4A in this case, so we do get an equation for the metric (note in 
pass the R =constant follows just from the scalar field equation (12bl) when = constant). 

In this paper we will consider only the case F^ u = 0, and will explore the space of static, 
spherically symmetric solutions: 

ds 2 = -N 2 (r) dt 2 + Nt(r) dr 2 + r 2 dtt 2 , = 0(r). (10) 

Since all solutions of the field equations satisfy PJ, we will oftentimes replace eq. ( l2bl) with 
the much simpler equation 

2 

□0 - -A0 - 4a0 3 = 0. (11) 
o 
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We were naturally led to consider this problem from the linear stability analysis in 20| of 
the exact solutions found in These exact solutions are all of the form 

ds 2 = -N{r)dt 2 + Nir^dr 2 + rW, = <f>(r). (12) 

The first one that we analyze, which we call here solution MTZ1, has a constant scalar field 
(JHJ). For generic theories N(r) is obtained by imposing on (|T2l the condition G^ u = —Ag^ u , 
which follows from (l5aj) , or more directly from (JTj) . For the special theories © , as explained 
above, the Euler-Lagrange equation for the metric is trivial, and the only constraint on (Tl2]) 
is R = 4A, and comes from the Euler-Lagrange equation (15bl) for the scalar field. Since this 
condition on the metric is less restrictive than the one for generic theories, we get a wider 
set of solutions for special theories (two integration constants, Q and M below, instead of 
one): 

N(r) = 1- — - ^r 2 , fa) = M,a^ -A/36 (13) 

To V 



2M Q A „ r- 
N(r) = 1 - — + - -r 2 , <p{r) = y/E, a = -A/36. (14) 

T o 



In other words, requiring R = 4A to the metric (|12p gives N(r) as in (|14p . Adding the extra 
condition R^ v = Ag^ forces Q = 0. Note that ([TBI is the Schwarzchild-(A)dS metric in the 
generic case, Reissner-Nordstrom (A)dS for special theories (with Q the "source" of a scalar 
field instead of the square of the electric charge). 

Since we are only interested in static black hole solutions with a physically acceptable stress- 
energy-momentum tensor, we require that the singularity at r = be hidden behind an event 
horizon, and that Tj% v satisfies appropriate energy conditions in the N(r) > region between 
the event and cosmological horizons, which we assume located at r > 0, using if necessary 
the invariance of the metric under (r, M) -> (-r, —M)). MTZ1 has Tff v = + Ag^ = 
in the generic case. A straightforward calculation shows that for the special theory ([TJ 



T H"=% ( V^' - V I> V " +Ulfi<> - °l>°!> ) <b r )) 
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in the natural orthonormal basis t M = N^ 1 ^ 2 dt, f M = N 1 / 2 d r , 6* M = r~ 1 dg 
(r sin(6')) _1 9 9 i. Thus, the strong and dominant energy conditions are satisfied in both cases 
as long as Q > 0. 

The second type of F^ v = solution in [8| for the system (pQ) and the Ansatz (TT2"]) . which 
we call MTZ2, holds only for the special theories a = —A/36. The metric is that of a 
Reissner-Nordstrom (A)dS black hole, the mass being an integration constant that appears 
both in the metric and the scalar field: 

^ ~ r-M 

To avoid naked singularities, we restrict to the case A > 0, then — > — ooasr — > oo, 
A^ — > oo as r — > + , and the singularity at r = is not naked only if N(r) has three positive 
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roots. This can only happen if < M < v / 3/(4\/A) =: 1/4. In this case, the three positive 
roots are 

/ / A 1\/T \ 1 I I A 1\/T\ 1 I I A 1\/T\ 

(17) 

This solutions are black holes on a cosmological background, with an inner horizon r\, a 
regular event horizon r 2 and a cosmological horizon r% [8|. 

It will be useful for our discussion to review the derivation of the MTZ metrics from the 
Ansatz f|T2|) . Notice that the Einstein plus scalar field equations imply that the functions 
N(r) and <p(r) must satisfy four equations, and therefore the set of solutions is severely 
restricted. Assuming as stated that F^ v = 0, and the form f[T2"l) for the metric, an appropriate 
combination of the Einstein equations implies that 4> satisfies the equation, 

• rt -2^V = (18) 



dr 2 \ dr r 

This admits the solution <j>(r) = 0, leading to vacuum black holes with A ^ 0, and also a 
general solution of the form, 

« r > = c-TTci (19) 

where C\, and C 2 are constants. The two kinds of solutions, MTZ1 and MTZ2, are obtained 
by choosing C\ = or C\ ^ 0, then solving the remaining field equations. There is no other 
solution to the field equations of the form flty) . 



As far as we know, a linear stability analysis of MTZ1 has not yet been done. In what 
follows we sketch the construction of some particular unstable modes for the theory (jSj), of 
the restricted form, 

S(f)(r,t) = 
Sg rr (r,t) = F(r) exp(kt) 

5g tt (r,t) = -AF(r)exp(kt) (20) 

where 5 indicates the perturbed part, and A is constant. Unstable modes would result if we 
find appropriate solutions for the perturbation equations with k real and positive. Replacing 
this Ansatz in Einstein's and the scalar field equation, and keeping only linear terms in F, 
the only non trivial equation that results is of the form, 

tl = M± d _l + F{r) + k> 

dr 2 r 7 N{r) dr r w N{r) 2 A 1 ' 

where Pi and P2 are polynomials in r with coefficients that depend only on A, M, Q, and A, 
which are therefore regular in the relevant range in r, that is for r# < r < r&, with r = r# 
(the event horizon), and r = r& (the cosmological horizon) corresponding to single zeros of 
N. It can be checked that, the general solution of (|21|) . near one of the zeros of N, which 
are the singular points of f|2T|) . is of the form, 



ci(r - r p ) + c 2 \ \r - r p \ (22) 
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where c\, and C2 are arbitrary constants, and r p is either r# or r\. This means that the 
general solution of (12TT) vanishes at the horizons, but only those solutions with C2 = at both 
r = th and r = are acceptable as perturbations, because for C2 7^ the derivatives of F(r) 
are singular. This implies that appropriate solutions, if they exist, satisfy a boundary value 
problem, with k 2 the corresponding eigenvalue. Considering now a numerical integration of 
(12T1) . there is no difficulty in imposing regularity for, say, r = r#, but, for general A and 
k the resulting solution would be singular for r = r\. We notice, however, that for large 
enough k, and r not close to the horizons, (12T]) behaves approximately as, 



d 2 F k 2 

Therefore, if we take A < 0, F(r) will oscillate between positive and negative values in the 
region r# < r < ta. This implies that, for A < 0, imposing the condition that C2 vanishes 
for both r = th and r = r\ turns (12TT) into a boundary value problem determining the 
allowed values of k. Note that (1221) guarantees that the perturbation will vanish at both 
horizons. Clearly, there is no upper bound on the allowed k values. Therefore, the linear 
perturbation problem leads to solutions that diverge arbitrarily fast from MTZ1. Figured] 
illustrates a "shooting" approach to the problem of finding appropriate values for k: Q,M 
and A were chosen so that r# = 2 and = 16, and (12"T|) was numerically integrated from 
r = th, setting C2 = at this horizon (see eq. ff22l ). Generically, the solution at r = r\ will 
also be of the form (1221) . but with C2 7^ 0, then F' will diverge there. Requiring that F' be 
finite at both horizons gives a discrete set of possible k values. The left panel of the figure 
shows a numerical integration performed with k = 1.0, the right panel shows a numerical 
integration with k = 1.2. The fact that at both horizons the behavior is as in (T221 guarantees 
the vanishing of F. It is clear from Figure [Hand continuity arguments, that, for some value 
of k in this interval there is a solution with a finite derivative at r^. 

The analysis carried out in [20j indicates that the solutions MTZ2 are also unstable 
under linear, spherical perturbations. Once again, if one attempts to solve the linear 
perturbation equations for the spherically symmetric mode, one finds solutions that grow in 
time arbitrarily fast [loj. This may be traced to the fact that the perturbation "potential" 
(of the Regge- Wheeler like equation) is singular for r = 2M, a rather peculiar situation, 
since the metric (Tl6l) and the scalar field are smooth in the range between the event and 
cosmological horizons, in particular at r = 2M, since r 2 < 2M < r 3 . 

More generally, the problem of solving the linearized equations for arbitrary (i.e., not 
restricted to spherically symmetric) perturbations can be approached by decomposing 
in angular modes in the usual way, and projecting onto 5* 2 harmonic vector and scalar 
fields, but, as we have checked, this leads to and extremely intricate set of equations that 
is difficult to deal with. Notice however that in order to prove instability, it is certainly 
sufficient to exhibit a single unstable mode, as was done above for MTZ1 and in [20j for 
MTZ2. 



An important point is that under the radial perturbations above and in 20], the 
perturbed metrics leave the restricted family g tt = —l/g rr , getting into the general space of 
static and spherically symmetric spacetimes ffTO]) . This suggests that the peculiar behavior 
of the MTZ solutions under perturbations may be related to the restricted nature of the 
space of solutions of the form (1T21) . For instance, the perturbation method might not be 
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FIG. 1: Numerical integration of equation (|2ip from r = r# to r = r\. Equation (|22p guarantees 
that F will vanish at both horizons, however, F' will generically diverge at r = r\ if it is finite at 
r = rn, except for a discrete set of k values. The allowed k values can be spotted by a numerical 
"shooting" algorithm. This is illustrated in the figure above: the left panel shows a numerical 
integration for k = 1.0, the right panel a numerical integration for k = 1.2. It follows that there 
is an allowed value k Q , with 1.0 < k Q < 1.2. The integrations were performed for a special theory 
([8]), setting A = Q = ^A, and M = which gives r H = 2 and r A = 16. 



applicable because in the general case N 2 ^ l/Aq, there are solutions which are locally 
arbitrarily close the unperturbed one in the family (ITUl) . but with very different global be- 
haviour. This, as we show in this Paper, is precisely the case for the MTZ family of solutions. 



III. THE EINSTEIN EQUATIONS 



We generalize the metric Ansatz of [8| by considering instead of (IT21 . a static spherically 
symmetric metric and scalar field of the form (TTUj) : 

ds 2 = -N 2 (r)dt 2 + Aq(r)dr 2 + r 2 dVL 2 , = <f>{r) 

Recall from the previous Section that (jl]) always hold, and thus ( I2bl) can be replaced with 
( TTTT) . Inserting ( TTOl) in (Aq times) Eq. ( TTTT) gives: 







V 



N 2 ' Ni 4' 
iV^ ~ iV\ + r 



-N x 4> [6a<p 2 + A] . 



(24) 



Also, Q^v := G^ v + Ag^ — T^ v is diagonal, with Qgg oc , thus ( |2al) gives three non trivial 
equations. The first two are: 







6r 2 N\ 
No 



AT/ 



- 6) + r<t>'<j>] - Nx [((j) 2 - 6) + 6r 2 (A + a0 4 )] 

+ [4r00' + (0 2 - 6) + 2r 2 00" - r 2 0' 2 ] , (25) 
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N ' 

= 6r 2 Q„ = --f r [(0 2 - 6) + r<f>'<f>] + N x [(0 2 - 6) + 6r 2 (A + a0 4 )] 

- [3r 2 /2 + 4r<f><f/ + (0 2 - 6)] , (26) 

The Goo equation is seen (after some work) to actually follow from fl24|) . fl25|) and fl26|) . so 
it will not be needed. The field equations conform a system of three ODEs, ((23]) -([26]), on 
three unknown functions N\, N2 and 0. 

It is apparent from (I2l)|) - fl2"6"|) that the case where [{(f) 2 — 6) + r0'0] = is special. If such 
a solution exists then = V6r 2 + C/r and equations (|25|) - (|26|) force C = (i.e., = v^), 
and a = —A/36. Under these conditions the remaining field equation, eq. (12^|) . is also 
satisfied. This is of course solution MTZ1 in the special case a = —A/36, eq. (TH|) . 
If, on the other hand, [(0 2 — 6) + r0'0] ^ (in particular 2 ^ 6), we find, after some work 
on (E4D-d26D that 



[(2 (9 A - <j) 2 A + 3 a 4 ) r 2 + 3 <f) 2 - 18) 0' - 2 r0 (0 2 - 6) (6 a 2 + A)] JVi 

£2) 

(27) 



<ir 2 3r (6 

(6-0 2 + (0')V + 2r00') 0' 



r (6 — 2 ) 

and also that we can write Tff u just in terms of N\, N%, and 0', using the orthonormal basis 
t> = A" 2 ~ 1/2 d t , r>* = iVr 1/2 d r , = r^de, & = (r sin^))" 1 ^: 

Tya, = p ij, v + p r f^f v + p e 0^ + p^ 0^. (28) 

Here 

(-12 5 r 2 a + (A r 2 - 3) 3 + (18 - 18 A r 2 ) 0) 0' + 6 a 6 r + (-36 ra + A r) 4 - 6 A r0 2 
9 ~ 3r (-6 + 2 ) (-6 + 2 + r0 0') 

6 r 2 0' 3 + (18 r + 9 r0 2 ) 0' 2 + (3 3 - 18 0) 0' 



Pr 



4 3r(-6 + 2 ) (-6 + 2 + r00') A; 

(-l+Ar 2 )00' + r0 2 (6a0 2 + A) 0' (r0' + 



r(-6 + 2 + r00') r (-6 + 2 + r00') A; 

_0 2 (6a0 2 + A) 0'(r0' + 20) 
Pe ~ n ~ 3(0 2 - 6) r(0 2 - 6)iVi 

These formulae do not hold, of course, at those isolated points where (—6 + r00' + 2 ) = 
or 2 = 6. The fact that this system is singular at points where 2 = 6 is clearly related to 
the singular nature of the linearized equations for perturbations of MTZ2 (Fl6|) at r = 2M, 
where = y/E 0. This is so because the perturbed metric (eq.(9)) and scalar field (eq. 
(11)) in [2cj] are, in the static case, of the form (flQl) . In principle, whether a solution is also 
singular at such a point depends critically on the behavior of the numerators of equation 
(|27l) . The exact solutions (TIB!) represent cases where this singularity is canceled, but other 
possibilities should be expected. 

For the MTZ1 solution in the special case a = —A/36, eq. (j!4p . we cannot use (1281) . and 
the energy- momentum tensor in this case is given by (fT5l) . For ffl3|) . using (|28p we get the 
expected result T^ v = 0. 
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IV. SOLUTIONS WITH A REGULAR HORIZON 

In this Section we consider solutions of the field equations that: (i) contain a regular 
event horizon at r = ro, with tq > 0, and (ii) satisfy the weak and dominant energy 
conditions in some open neighborhood < r < r + e outside the horizon. 

The regular horizon condition implies that there exists a neighborhood of r = r , where 
the functions N\, N 2 and admit expansions of the form, 

= a + ai(r - r ) + a 2 (r - r ) 2 + . . . 
Nt = b-i(r - ro)- 1 + b + h(r - r ) + . . . (29) 
N 2 = ci(r - r ) + c 2 (r - r ) 2 + c 3 (r - r ) 3 + . . . 

where a», 6j, and Cj are constants coefficients. The proper signature of the metric imposes 
b-i > and ci > 0, although c\ is otherwise arbitrary because of the freedom of rescaling of 
t. We also impose ao > 0, making using of the invariance of the equations under — > —0. 
We will find it convenient to introduce the dimensionless horizon radius 

z := r VX. (30) 

In Section IIVBI we arrive at a description of the subset of the (z a , a ) plane for which 
conditions (i) and (ii) above are satisfied in the special case a = —A/36 (see Figure |2j). 

Replacing the expansions ( 1291 in equations (I24l) -( l26l we obtain relations between the 
coefficients by equating powers in r — r . From the algebraic equations obtained by 
matching the lowest order non trivial terms we learn that the "special" case a = —A/36 
requires separate treatment. 



A. Generic theories (a / —A/36) 
To lowest order we obtain 

6-i 



2r a (a 2 — 6) (6 a a 2 + A) 



3 a 2 - 2 a 2 A r 2 + 6 a a 4 r 2 - 18 + 18 A r 2 ' 
-3r (-6 + a 2 ) 



18-3 a 2 - 18 A r 2 - 6 a a W + 2 a 2 r 2 A 
This suggest that we study the cases A) a\ — and B) a% ^ separately. 



(31) 
(32) 



Case A.l, clq = 0: in this case, by solving iteratively for the higher order terms, assuming 
r > 0, we are led to the Taylor expansion of Schwarzschild - dS (S-dS) space: 



(r) = 0, Ni(r) 



1 - 



r (3 - Ar 2 ) Ar' 



3r 



N 2 (r) 



CiTq 



(Ar 2 - l)iVi(r) 



(33) 



Matching N\ — 1 — 2M/r — Ar 2 /3 = — ^(r + r x )(r + r 2 )(r + r 3 ), gives r x = — (r 2 + r 3 ). We 
want, say r 2 = ro (event horizon), r 3 = cosmological horizon, then < r 2 < r 3 . This implies 
M = A(r 2 + r 3 )r 2 r 3 /6 > and 

z < 1. (34) 
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Case A. 2, a = Vd: this leads to b_i = 0, see eq. (1521) . Interestingly, no solution with a 
regular horizon and 0(ro) = Vd exists in the generic theory if we assume the scalar field 
admits a Taylor expansion around the horizon. 



Case A. 3, a = w-^: in this case we obtain, once again, aj = for j > 0, i.e., 
together with 



A 
6a ' 



NUr) 



r 



1 - Ary 



+ 



1 + (1 - Ar 2 
Ar (4 - Ar 2 



r ) + 



r (3Ar 2 + 1 - A 2 r 4 



3(1 - Ar 2 ) 3 v ' 3(1 - Ar 2 ) 4 

which is the Taylor expansion around r = r of MTZ1, eq. ffTB"]) . written as 

r (Ar 2 - 3) Ar 2 



-(r - r ) 2 + . . . 



1 + 



3r 



-i 



Case 5, ai 7^ 0: This case is extremely complex to deal with in the general situation. 
Since the main motivation of this work is to understand the behavior of the linearized field 
equations around the solution MTZ2 ffT6l) for the theory a = — ^, we restrict our attention 
to special theories from now on. 



B. Special Theories a = —A/36 

To lowest order, for the special theories we obtain 

-2Aroao(a 2 - 6) = 18ro 

fll ~ 18-r 2 A(a 2 + 18)' ^ ~ 18 - r 2 A(a 2 + 18) ' 1 ' 

This suggests that we study the cases ao = and a = VE separately. 

Case A: a Q = 0: To no surprise, we are led back to Schwarzschild de Sitter space, eqs. f )33|) 
and (El. 



Case B: a = The higher order terms of ( |24|) . ( 1251) and ( 1261) imply give a-,- = 0, j > 0, 
i.e., any solution with 0(ro) = \/6 must satisfy 0(r) = \/6 for all r. From the comments in 
Section [III we know that the only field equation for the metric in this case is R = 4A, which 
reads 

n 2 " (iW) 2 aw 2ivy 2ivy 2^-1) _ 



A^Aa 2NiN 2 2JVJV2 rNtN 2 rN x 2 r 2 ^ 

In principle, this gives us an infinite number of solutions for the Ansatz (TTU|) . since, given, 
say A^ 2 , R = 4A is a first order ODE for N\. In particular, given A" 2 as in f)29p and any 
6_i > 0, the algebraic equations for the remaining admit a solution. Inserting this 
solution in the energy momentum tensor T^ v = G^ u + Ag^ v and using the orthonormal basis 
in fT28|) gives 

, <ntt \\ 6-i(l -Ar 2 ) -r 
p = p o + C((r-r )), p Q := — — , (37) 
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and 



Pe Pr , . /2ci6_i(2Ar 2 -l)+4ciro + 2c 2 ro 2 \ , , , 2 x , , 

— — = 1 + F ; — 71 — o Tvi (r - r ) + C (r - r ) (38) 

The condition p D > is equivalent to 

l-Zo 2 1 _ y/l+AAb^-1 

> — == — <^> < z Q < == , (39) 

Zo VAb-i 2VA&-1 

and thus z < 1, as happens for generic theories, eq. (134]) . To satisfy the strong energy 
condition in some open r interval r < r < r + e we require that the (r — r ) coefficient 
in (13"8"|) be negative, and this can always be satisfied by a proper choice of c 2 , thus proving 
that there are local solutions satisfying the energy conditions right outside a regular horizon. 

It is not hard to see that, out of the infinitely many solutions for the ODE fl36l) . the 
only one satisfying N\N 2 = 1 is MTZ1 f fl"4l) . Given MTZ1 with positive A and positive 
Q (required by the energy conditions, see ffl5l) ). one can easily see that in order to avoid 
naked singularities the quartic polynomial r 2 N(r) has to have four real roots, one negative 
and three positive: — < < 7*1 < r 2 < r 3 , with r 2 = r$ the event horizon and r 3 the 
cosmological horizon. Then matching (|14p with 

N(r) = --^( r + r i)( r - r i)( r - r 2)(r - r 3 ), (40) 

gives r 4 = (ri + r 2 + r 3 ), a positive mass 

M = (r 1 +r 2 )(r 1 + r 3 )( y r 2 + r 3 ) ^ 

and 

\ qr> qr> \ ry ry 

l^i<j<Z'v] A~ti<j<3 1 3 

In particular, given the domain < r x < r 2 < r 3 < 00, one finds that 2 = r 2 \/A satisfies 
the constraint (compare to (jM]) and (157)) ) 

< z < 1 (43) 

Case C: ao 7^ 0, Vd: In this case we find that all the coefficients in (1291) may be written, 
e.g., in terms of r , and a . The leading terms are of the form, 

2(a 2 - 6)Ar a 4(a 2 - 6) 2 A 2 r 2 a , , 2 

" 00 + Ar 2 a 2 + 18Ar 2 - 18 ( " " ro) + (Ar 2 a 2 + 18Ar 2 - 18) 2 ( " " r ° } + '- 

Nl = ^ 36 ( A « ' g) + (44 ) 

1 (18- 18Ar 2 -Ar 2 a 2 )(r-r ) (18 - 18Ar 2 - Ar 2 a 2 ) 2 1 } 



N 2 = Cl 



^'o^o 

2„2 



, . 2(Ar 2 a 2 -9) , , 2 

(r -''° ) + M18-18A°g-Argag) (r - r ° ) 



where, as already noticed, C\ > 0, but it is otherwise arbitrary. This implies that the 
condition for the existence of a regular horizon leads, in general, to a two-parameter (ro 
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and do) family of solutions. We notice, for reference, that the exact solution MTZ2 (Fl6|) 
corresponds to the one parameter subfamily for which 

r = — f= : 1 o = , < M < - W — (45) 

2v/A 2A ' r -M' - -4VA 



with ci chosen as Ci = ( ' 18 18 ^g° o Ar ° a °' > in (ITBl) . The limit case M = gives just de Sitter 
spacetime with no scalar field. 

As explained in Section [Til MTZ1 and MTZ2 are the only solutions with N% = N 2 . It 
is important to check that the expansions (j44l) are consistent with this fact. From (|44p we 
obtain 

N N = 18r o c i 24ciA7Wg _ 2 

1 2 ~ 18 - r 2 (a 2 + 18) A [_1 8 + r 2 (a 2 + 18) A] 4 1 ° J 
_16ciAa 2 [Ar 2 (a 2 - 66) + 9] B , 4 

[-18 + r 2 (a 2 + 18)A] 5 1 ° J U °^ 1 J 

where, 

B = 324 + r 4 (a 2 - 6) 2 A 2 - 36 r 2 (6 + a 2 ) A (47) 

Therefore, the condition Ni(r)N 2 (r) = 1 can be imposed only if a = 0, which is trivial, or 
if B = 0. In this case, solving for a 2 , in terms of the other constants, we find two solutions, 
but only one of these leads to acceptable coefficients in (1441) . This solution is given by, 

2 18 + 6Ar 2 -12V3r VA 6(^3 - z f 

a = TZ2 = — 2 I 48 ) 

11 I z o 

and it can be checked that this coincides with (1431). 



Another interesting issue is that of analyzing the limit a — > v6 in (T44|) . The limit gives 
= v^6, and well defined expansions for iVi, and jV 2 , that can be seen to satisfy the required 
condition on the metric, R = 4A. Thus, this is one of the infinitely many = solutions 
referred to in Case B above, certainly not MTZ1 (TT4"]) . since do — > \/6 in (g6D gives 

iVx(r)iV 2 (r) = 3 gg^2 + (^-^ ( r " ^ + ° (( r " r °) 3 ) + constant ( 49 ) 

We may obtain important information regarding the physical acceptability of the solutions 
( T44|) by considering the behavior of the energy-momentum tensor near the horizon. Imposing 
the strong and dominant energy conditions on (|44j) places restrictions on the range of the 
parameters (r , ao)- I n the notation of equation 



1 ?\n 2 

P=^Aa 2 -^(r-r ) + O((r-r ) 2 ), (50) 
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and 



Pr 
P 



£Z = _i 



Ar 2 a 2 - 6( v / 3 + rov/A) 2 Ar 2 a 2 - 6( v / 3 - r VXf 



3r 2 (Ar 2 a 2 - 18 + 18Ar 2 ) 2 



(r - ro) : 



+C((r-r ) 



4 



Ar 2 a 2 - 6(V3 + r VA) 2 Ar 2 a 2 - 6(^3 - r x/A) 5 



1 + 



+0((r-r o ) 3 ) 



3r 2 (Ar 2 ag - 18 + 18Ar 2 



2\2 



(r - r f 



(51) 



Therefore, the solutions satisfy the weak energy condition, (positive energy density) for 
all ao, but they violate the dominant energy condition (absolute value of the stresses not 
larger than energy density) in the neighborhood of the horizon unless ao and ro are restricted 
by the conditions, 

6(v / 3-r v / A) 2 < Ar 2 a 2 < 6(v / 3 + r v / A) 2 (52) 

At the limits we have p r /p = — 1, and Pe/p = 1- The upper limit is further restricted by the 
condition, 

Ar 2 a 2 < 18 - 18r 2 A (53) 
imposed by the condition Ni > 0. All together this implies, 



6(^3 -z ) 2 2 18(1 -z 2 

o S a < 5 



(54) 



Note from (145]) . (1481) . that MTZ2 ([16]) saturates the lower bound above, and that the allowed 
interval for a 2 is nonempty only if 

z <^y- (55) 
The restrictions for case C can then be summarized by any of the two equivalent conditions: 



6 < 



6(v / 3 



< ao < 



18(1 - ^ 



(the first bound in the chain of inequalities following from (1351) ). or 



1 4- -22- 



<z < 



< 



(56) 



(57) 



+ 3 



(the last bound in the chain of inequalities following from ao 2 > 6). This completes the 
discussion of case C. 



Let us recapitulate on what we have found by seeking local solutions of the form (129]) 
for the special theories a = —A/36, satisfying the weak and dominant energy conditions 
outside the horizon. We have used the <fi — > —0 symmetry of the field equations to restrict 
our considerations to 0(r o ) =: a > and found that: 

• If ao = then <p(r) = and the metric is Schwarzschild de Sitter. The constraint 
z < 1 is required to assure there is an event horizon hiding the singularity, and an 
exterior cosmological horizon. 
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• For < ao < v6 there are no solutions satisfying the weak and dominant energy 
conditions outside the horizon. 

• If a = then <p{r) = v^6, and there are infinitely many solutions (Ni(r), N 2 (r)), for 
every z D < 1, (see eq. ([5rJj) ). some of them satisfying the energy conditions. 

• If ao > v^6 then there is one solution satisfying the desired energy conditions for every 
pair (ao, z Q ) satisfying ( JB"7[) 

The situation is summarized in Fig[2j 




FIG. 2: Allowed regions in the (z ,a ) plane for static spherically symmetric local solutions of the 
special theory a = —A/36 having a regular horizon, eqs. (jlOp (|29p and satisfying the weak and 
dominant energy condition in some exterior neighborhood of the horizon. There is one solution 
per point except at the subset clq = V6, < z < 1, where there are (infinitely) many solutions 
per point. The horizon radius is ro = \/A z q and the value of the scalar field at the horizon is 
0(ro) = ao- The ao = 0,z o < 1 solutions are Schwarzschild - de Sitter, the MTZ1 solutions lie on 
the segment ao = V6, < z < 1, the MTZ2 solutions on the lower edge z Q = of the shaded 

VD+ao 

ao > V® region. 

The natural question to ask at this point is what is the global behavior of the local solu- 
tions analyzed above. Since solving the system f|24|) - f|26]) analytically is out of consideration, 
numerical integrations were performed. The results are gathered in the following Section. 

V. SPECIAL THEORIES: NUMERICAL ANALYSIS OF THE 0(r o ) > V6 

SOLUTIONS 

The equivalent conditions given in equations ( 1561) and ( l5"Tj) provide a range of values 
for r and a such that, locally, the field equations have a solution with a regular event 



16 



horizon, with the strong and dominant energy conditions being satisfied right outside the 
horizon. The question that naturally arises then is what is the behavior of these solutions 
as we move away from r = r . Since we do not know of any exact solutions in this range 
besides the borderline MTZ2, we considered a numerical integration of the system (1241) -( 1261) , 
using the expansions ( 1441) to construct appropriate initial data. A numerical integration 
requires assigning definite numerical values to the parameters. We took C\ — 1, A = 3, 
and r = 1/4, and considered different values of ao in the interval 3v^6 < ao < V78 (eq. 
(156]) ). To check the accuracy of the numerical procedure [24]], we analyzed as a first example 
the MTZ2 data ao = 3\^Q, which corresponds to M = 3/16 in eq. (fTBI) . With this choice 
of ao we have <fi 2 = 6 for r = 3/8, and the equations are formally singular, because of 
vanishing denominators, for this value of r, with the result that the numerical integration 
stops at that point. Nevertheless, the numerical solution is well behaved for any r close to 
but smaller than 3/8, in correspondence with the regularity of the exact solution, with a five 
digit agreement between the exact and numerical solutions in the plotted range, Figure [33 




FIG. 3: Numerically generated MTZ2 solution with A = 3,ro = 1/4, ao = 3\/6. The vertical axis 
displays the correct values of the scalar field, the scales of N\ and N2 are arbitrary, and were 
independently chosen to fit the range of 4> values. Notice the smooth behavior as <j) approaches the 
regular singularity at <fi = \/6. 

Next we considered, for the same value of r , a number of different allowed values of ao 
larger than 3\^6. The general behavior turned out to be qualitatively the same in all cases: 
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the numerical integration shows a singular behavior in N\, as r approaches a critical value 
r = r$, while A^ 2 and <p approach finite limits, with <p - ► 0c > v6- This is illustrated in 
Figure HI It is also found numerically that T t f , as well as other invariants, approach a finite 
limit as N% diverges. This raises the possibility that the singular behaviour for r = rs is 
only a coordinate effect. In the next Section we show that this is effectively the case. 




FIG. 4: Numerically generated solution with A = 3,ro = 1/4, do = 3\/6 + 0.3. The vertical axis 
displays the correct values of the scalar field, while the scales for N\ and N 2 are arbitrary, and 
were independently chosen to fit the range of <j) values. As r — > 7*5, approaches a critical value 
4> c > V® and N\ diverges. The critical value 4> c ~ y/E in this example because ro, and ao are close 
to the MTZ2 values. 



A. A coordinate singularity and extensions of the solutions 

A detailed numerical analysis of the behavior of N%, N 2 , and <fi near the singular point 
r = rs indicates that, in general, we have, for r < rs, and r ~ rs 

iVi(r) Ui 



rs — r 



N 2 {r) ~ + 

~ A + AiVrs-r (58) 
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where, Ao, Ax, Bx, Co, C\ are constants that depend on the solution, and A > Vo. This 
suggests the introduction of a new coordinate R, defined by, 

R = v^^7 (59) 

For this new coordinate the metric takes the form, 

ds 2 = -N 2 dt 2 + NidR 2 + (r 5 - R 2 ) 2 dQ 2 (60) 

where 

iVi = 4R 2 iVi 

N 2 = N 2 . (61) 

The resulting Einstein and scalar field equations in the new coordinate R are rather long, 
and we do not display them here. We find that, just as in the case of the r coordinate, 
they are equivalent to a set of three equations for Nx, N 2 , and 0. The system has singular 
coefficients for R = 0, but admits regular solutions in the neighborhood of R = 0, with Nx, 
N 2 , and <ft having expansions of the form, 

Nx = B + B X R + B 2 R 2 + ... 
N 2 = C + CxR + C 2 R 2 + ... 
<f) = A + A 1 R + A 2 R 2 + ... (62) 

where B^, and Ci are constants, and B = 4Bi,C = C ,Ci = C\, in agreement with 
( 158]) . Since the transformation fl59|) is defined only for R > while (1601 is defined also for 
R < 0, the coordinate change (1591) provides a smooth extension of the original metric (Tip"]) 
through the singular point r = rs- 

We are again here confronted with the lack of explicit exact solutions, and, therefore, 
we must resort to a numerical integration to obtain information on the properties of these 
solutions. This is considered in the next Section. 



B. Numerical analysis of the continued metrics 

Given the form (160]) for the metric, a regular horizon at R = Rh = ^/r^ — r would be 
characterized by the functions N\, N 2 , and <fi admitting expansions 

h , 



Rh — R 

N 2 = cx(R H -R) + d 2 {R H -R) 2 + ... 
= a + ax(R H - R) + ■■■ (63) 

where Sj, bi, and Cj are constants. Given a particular solution for (TlOl . with a regular horizon 
characterized by given values of ro, and ao, and the remaining coefficients given by ( |44"j) . for 
which the singularity appears at r = rs, we have the following relations for the coefficients 
of the leading terms: 



Ci = 2c lv / r s - r (64) 
36r 0N /r s - r 
~ 18- (a 2 + 18)r 2 A 
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We use as initial data for a numerical integration of the equations for iV 1; N 2 , and 
in the region Rh > R > 0. The numerical integration stops for R = 0, but shows that 
iVi, A^2, and <fi display a regular behavior arbitrarily close to R = 0, and allows to extract 
the leading coefficients Ai, Bi, and Cj in ( 1621 to compute initial data for the numerical 
integration of the equations in the region R < 0, i.e., beyond rs- 

The main result is that in this extension we find that <p — > y/E, while N% — > 0, and 
A^ 2 — > oo as i? — > — ^/rg. The energy density p, is found to diverge as i? — > — \ffs- 

This situation may be analyzed in general by noticing that R — > —y/rs corresponds to 
r —>■ + in pop if we change variables to r = — R 2 , so that r — > as i? — > — v/rg. The 
numerical results suggest that (0 2 — 6) — > 0, and N% — > as some power of r, while both 
A^2 and p diverge. The detailed behavior near the singularity depends, however, on some 
rather delicate cancelations of diverging terms, and, up to the accuracy achieved so far, 
we can only draw qualitative conclusions out of the numerical results [25|. To this extent, 
it appears that the extensions end at a (naked) singularity, and that the solutions cannot 
be further extended. This would imply that the only solution with ao > v6 containing a 
region limited by event and cosmological horizons, where the energy momentum tensor is 
compatible with the weak and dominant energy conditions, is the exact solution MTZ2 
found in [§]. 

Nevertheless, for the problem of understanding the instability found in [20j we need to 
study the neighborhood of the MTZ2 curve in Fig [21 and this includes the dominant energy 
violating cases where 6(\/3 — rov^A) 2 > Ar^a 2 ,. These are considered in the next Section. 

C. Solutions violating the dominant energy condition 

A numerical analysis of solutions with 6(v / 3 — z ) 2 /z 2 > a 2 , -i.e., violating the dominant 
energy condition near the horizon and thus (156jl - reveals a smooth behavior of the metric 
for r > r , with (j) — ► Vd, and N\ — > + , N 2 — > + as r increases past some value larger than 
tq. The most remarkable feature of these solutions is that the energy density p(r) decreases 
from its value at the horizon r = Tq, changing sign at some r\ > ro, with p taking larger 
and larger negative values as r increases. As already mentioned, the numerical integration 
breaks down for sufficiently large values of r, but not before the divergence of \p\ is clearly 
established, leading to the conclusion that solutions outside the allowed regions shown in 
figure [2] contain features that make them physically unacceptable. 

It is interesting that when a is slightly smaller than the lower bound forced in 
which, as we said, corresponds to the MTZ2 solution (fl6l) . the numerically generated solution 
remains close to the MTZ2 solution for r ~ ro, and then they depart completely from each 
other as we move away from the horizon. This is a coordinate independent statement, since 
it is exhibited, e.g., by a qualitatively different behavior of the energy density in both cases. 

VI. SUMMARY AND CONCLUSIONS 

We have studied the theory ((T]) with = 0, A > 0, and a = —A/36, and arrived 
at a comprehensive understanding of the space A4 of static, spherically symmetric local 
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solutions with a regular horizon that satisfy the strong and dominant energy conditions 
in an open set bounded by the horizon. The diagram in Figure [2] shows the (ao,z ) plane 
(z Q := a/Ato, ao = 0(po), 7*o the horizon radius.) We have proved that there is a one to one, 
onto correspondence between the set of ao 7^ a/6 solutions in Ai, and the (do 7^ a/6) shaded 
region of this plane. Among these, the only known exact solutions are MTZ2, eq. ffl6l) and 
= Schwarzschild - de Sitter spacetime. The case a = v6 is rather peculiar, for every 
point in the segment a = v6, < z Q < 1 there is not just one, but infinitely many local 
solutions of the field equations admitting a regular horizon, some of them satisfying the 
weak and dominant energy conditions. To this set belongs the other known exact solutions, 
MTZ1 given in eq. (TT4l). 

Numerical integrations of the field equations away from the horizon, indicate that those 
solutions in the ao > a/6 shaded area are not physically relevant, since they develop a 
singularity with infinite energy density, not protected by a horizon. It is rather interesting 
that, between this singularity and the horizon, a coordinate singularity was numerically 
spotted, and appropriate new coordinates could be constructed to cross over it. The spheres 
of symmetry (i.e., the orbits of the SO (3) isometry group) have a radius (square root of 
(47r) _1 times their area) that grows from the horizon radius r$ up to a maximum value r$ 
(where the coordinate change is required), and then collapses to zero as we approach the 
above mentioned (spacelike) naked singularity. 

The unshaded lower region in the a > a/6 portion of the (a , z Q ) plane corresponds 
to uninteresting solutions of the field equations. They not only violate the strong energy 
condition near the horizon, but also have an energy density p that, as we move away 
the horizon, becomes negative, and apparently unbounded as r increases, (no coordinate 
change is needed for these solutions). 

One of the main purposes of the present work was to obtain an understanding for the 
extreme instability under perturbations found in [20[ for the metric (1161) . From a simple 
perspective, given the family of solutions (|T6l) . one would expect that under a sufficiently 
small perturbation the system would radiate some gravitational and scalar field energy, 
partly to each horizon, and eventually settle to a static solution of the type (fT6|) . perhaps 
with different values of r and a (or AI in the notation of [8]). and therefore, the results 
of jlo| appear as difficult to understand. The present analysis, however, indicates that the 
parameter space for the static spherically symmetric solutions of the MTZ system indeed 
presents a sharp discontinuity at the exact solution, with neighbouring solutions displaying 
properties that depart completely from those of the solution (Tl6[ K In particular, the analysis 
of Section l"V AI shows that the coordinate system used both in [8J and [20j is inadequate for 
the perturbative study, because of the coordinate singularity intrinsic to that system. But 
the same analysis shows that even if the coordinate singularity is avoided, there are solutions 
that approach arbitrarily close to (fT6l) near the black hole event horizon at r = r , but then 
depart from each other with totally different geometrical properties. In fact, in accordance 
with (H8"|) . for a given A, the MTZ2 solution is obtained only if ao, r , and A are "fine tuned" 
so that (J3HD is satisfied, and any departure from that relation leads either to solutions with 
a divergent behaviour for finite r (before a cosmological horizon is reached), or to solutions 
with no cosmological horizon, but with a divergent behaviour for the energy density. 

The final conclusion of our analysis is that there appear to be no physically acceptable 
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stable solutions of the MTZ system that can be interpreted as black holes with a cosmological 
horizon in the exterior of is event horizon. 
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